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Abstract. We consider the nonrelativistic four-boson system in two dimen- 
sions interacting via a short-range attractive potential. For a weakly attractive 
potential with one shallow two-body bound state with binding energy B2, the 
binding energies Bpj of shallow iV-body bound states are universal and thus do 
not depend on the details of the interaction potential. We compute the four- 
body binding energies in an effective quantum mechanics approach. There 
are exactly two bound states: the ground state with B^ = 197.3(l)i?2 and 
one excited state with B^ = 25.5(1)-E>2- We compare our results to recent 
predictions for TV-body bound states with large N S> 1. 

Introduction. — Low-dimensional systems are ubiquitous in many areas of physics. 
In condensed matter physics, e.g., they are of interest in connection with high- 
ly superconductivity and the Quantum Hall effect. A two-dimensional boson 
system has been realized experimentally with hydrogen adsorbed on a helium 
surface [Q. The experimental progress with ultracold atomic gases and Bose- 
Einstein Condensates has made it possible to engineer low-dimensional atomic 
systems in anisotropic traps [H |3J 0] . 

This experimental progress has also stimulated theoretical activities. Re- 
cently, the problem of ./V weakly attractive bosons in two spatial dimensions (2D) 
was revisited 0. Relying on the asymptotic freedom of nonrelativistic bosons in 
2D with attractive short-range interactions, some surprising universal properties 
of self-bound iV-boson droplets were predicted for N S> 1. In particular, the size 
Rn of the A-boson droplets was found to decrease exponentially with A: 



Rn 



0.3417, A»l, 



(1) 



Rn-i 



*On leave from FZ Jiilich, Institut fur Kernphysik (Theorie), D-52425 Jiilich and HISKP (The- 
orie), Universitat Bonn, Nufiallee 14-16, D-53115 Bonn, Germany 



2 



Universal Properties of the Four-Boson System in Two Dimensions 



while their ground state energies Bjy increase exponentially with N: 



Bn 



8.567 



iV>l 



(2) 



Bn-i 



which implies that the energy required to take out one particle from the droplet 
is about 88% of the total binding energy. This is in contrast to typical three- 
dimensional systems, such as 4 He droplets, where the one-particle separation 
energy is much less than the total binding energy jSJ. Corrections to Eqs. lfT|) 
and (J2J are expected to start at order 1/N 2 and thus should be small for large 
N. For realistic interactions with a finite range ro, the equations are valid for TV 
large, but below a critical value, 



At N = iVd-it the size of the droplet is comparable to the range of the potential 
and universality is lost. If there is a large separation between R2 and ro, then 
iV cr }t is much larger than one. The universal regime also breaks down when the 
binding energy of the iV-body states approaches the same order as the particle 
mass and the particles become relativistic. For typical atomic and molecular 
systems, however, the breakdown of the universal regime will be set by the finite 
range of the interaction as specified in Eq. (j3J) . 

In light of the experimental and theoretical interest in weakly bound iV-body 
clusters in 2D, it is worthwhile to calculate their binding energies for short-range 
interactions explicitly. For a sufficiently shallow two-body state (such that the 
zero-range approximation can be applied), the binding energies for the three- 
body system in 2D are known exactly. They were first calculated by Bruch and 
Tjon in 1979 [7j. More recently, they were recalculated with higher precision [S1[S|: 
the ground state has B^ = 16.522688(1) B2 and there is one excited state with 
Bjp = 1.2704091(1)^2. The value of the ground state energy B^ 0) /B 2 differs 
from the prediction in Eq. (J2J) by a factor of two, indicating that the three-body 
system is not in the asymptotic regime. Of course, such deviations from Eq. © 
are expected for small values of N. However, the exact iV-body ground state 
energies should approach the prediction (J3J) as N is increased. 

The Four-Boson System in 2D. — In this paper, we calculate the binding 
energies of the four-body system with attractive, short-range interactions in 2D. 
We choose an effective quantum mechanics approach to generate and renormalize 
an effective zero-range interaction potential that reproduces a given two-body 
bound state energy. This potential will then be used in the Yakubovsky equations 
to compute the binding energies of the four-boson system. For a sufficiently 
shallow two-body bound state, the four-body binding energies for our effective 
potential and any realistic finite range potential will be the same. In this sense, 
our results are universal. This method was recently applied to the four-boson 
system in three spatial dimensions • It exploits a separation of scales in physical 
systems and is ideally suited to calculate their universal properties. In principle, 
finite range corrections to the universal results can be calculated systematically 
within this approach but they are beyond the scope of this work. 



l«iV< iV cr it « 0.931 In — + 0(1) . 



(3) 



L. Platter, H. W. Hammer, Ulf-G. Meifiner 



3 



In two dimensions, any attractive pair potential has at least one two-body 
bound state. If this bound state is sufficiently shallow, B 2 <C 7l 2 /( mr o)> its prop- 
erties are universal and the "true" interaction potential can be replaced by a 
(5-function in position space. In momentum space, such a 5-function potential 
corresponds to the interaction 



V(p,p')=A 2 , A 2 <0. 



(4) 



which is independent of the relative momenta p and p' of the incoming and out- 
going pair. The bare coupling constant A 2 is negative for attractive interactions. 
For convenience, we work in units where the mass M of the bosons and Planck's 
constant h are set to unity: M = h = 1. 

The interaction (j3J) is separable and the Lippmann-Schwinger equation for 
the two-body problem can be solved analytically. The S-wave projected two-body 
t-matrix t(E) is given by 



t{E) = 2vr 1/A 2 - 2vr 



dqq 



E 



-1 



(5) 



where A is an ultraviolet cutoff used to regulate the logarithmically divergent 
integral. The bare coupling A 2 = A 2 (yl) depends on the cutoff A in such a way 
that all low-energy observables are independent of A. Demanding that the t- 
matrix © has a pole at the bound state energy E = —B 2 , we can trade the 
coupling constant A 2 for B 2 and obtain the renormalized t-matrix 



t(E) = 2 



log 



Eh 

' E 



(6) 



which will be used in the following. 

The four-body binding energies can be computed by solving the Yakubovsky 
equations ^U|, which are based on a generalization of the Faddeev equations 
for the three-body system. The full four-body wave function is first decomposed 
into Faddeev components, followed by a second decomposition into Yakubovsky 
components. In the case of identical bosons, one ends up with two Yakubovsky 
components ipA and ipB- We start from the Yakubovsky equations in the form 
given in 



xPa = G t 12 P 

i> B = G t 12 P 



(1 + P 3 4)ipA + i>B 
(1 + P u )lP A + lfj B 



(7) 



where Gq = [E — Ek in \ is the free four-particle propagator and ii 2 is the two- 
body t-matrix in the subsystem of particles 1 and 2. The operator P^ exchanges 
particles i and j, while P and P are given by 



P — P13P23 + Pl2P r , 



12-T23 1 



and 



P = PiP 



13-T24 



(8) 



Including only S-waves and following the same steps as in Ref. jJJ], we derive a 
momentum space representation of Eqs. (J7J) in two spatial dimensions. We obtain 
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a system of two coupled integral equations in two variables: 

f°° f 27T dd> r 

F A (u 2 ,u 3 ) = 2t(E - \u\ - \u\) / du' 2 u' 2 / — G (u 2 ,u' 2 ,u 3 ) F A {u' 2 ,u 3 ) 

Jo Jo Z7r L 

f 2n d<f>' ~L 

+ / — { G o(u 2 ,U 2 ,U 3 ) F A (u 2 ,U 3 ) + G (u 2 ,V 2 ,V 3 ) Fb(v 2 ,V 3 )}{9) 
Jo Z7T J 



F B (v 2 ,v 3 ) = 2t(E - \v\ - v\) I dv[v[ ±G (v 3 ,v 2 ,v[)F B (v 2 ,v' 1 ) 







+ j G (v 3 , u 2 , u 3 )F A {u 2 , u 3 )J , (10) 



where we have used the abbreviations 



G ( Ul ,u 2 ,u 3 ) = [E - u\ - \ul - lul] , 

2 1„,2 „,21 _1 



G (v u v 2 ,v 3 ) = [E - v\ - \v\ - ug] , (11) 



and 



M 2 = \l\ u 2 + ff U 3 + «2«3 COS0' , -U 3 = Jll^ + - ^U 2 U 3 COS 0' . 



£ 2 = \/u 2 + |tt| + 1 1*2^3 cos^' , z} 3 = + - |n 2 U3 cos 0' . 



«2 = v/fl>2 + I u| - |u 2 f 3 COS 0' , U 3 = Jjvj + u| + V 2 V 3 COS 0' . 



^2 = y 4M2 + ^2 + u 2 u' 2 COS . (12) 

The integral equations © and ()10|) are well behaved and can be solved numeri- 
cally without introducing a regulator. 

Results. — The four-body binding energies can be found be discretizing 
Eqs. ((HI irU|) and calculating the eigenvalues of the resulting matrix. The en- 
ergies are given by the energies E < for which the matrix has unit eigenval- 
ues while the wave function is given by the corresponding eigenvectors. Similar 
to the three-body system, we find exactly two bound states: the ground state 
with = 197.3(1) 5 2 and one excited state with = 25.5(1) B 2 . These 
binding energies are universal properties of any four-body system in two spatial 
dimensions with an attractive short-range interaction and a sufficiently shallow 
two-body bound state with binding energy B 2 . Corrections from finite range ef- 
fects and three- or four-body forces are expected to be suppressed by powers of 
TqE where E is the typical energy scale of the problem. Therefore, finite range 
corrections will be more important for the deeper states. 

The binding energies for the four-body system in 2D have previously been 
studied for finite-range potentials using the Yakubovsky equations ^21 as well 
as variational and Monte Carlo methods |1.3| I14j. In these calculations the two- 
body binding energy was varied by about a factor of ten. Over the whole range 
of binding energies, the ratio B^/B^ was found to be approximately 2.9, in 
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Figure 1. Bn/Bn-i as a function of TV. The dot-dashed line shows the asymptotic value of 
8.567. The dashed line is an estimate of how the large- TV value is approached. 



apparent disagreement with our universal result. However, one has to keep in 
mind that the results from these explicit finite range potentials are far from 
the universal limit. The two-body binding energy in these calculations was not 
small enough for the zero-range approximation to be applicable. From Table III 
of Ref. e.g., it is clear that the ratio B% /B% is still a factor two away 
from the universal zero-range result even for the smallest value of B% considered. 
Furthermore, Fig. 2 of the same reference shows that the dependence of £3 /B2 
on B2 becomes very rapid near B2 = 0. Because the four-body states are more 
deeply bound, it is natural to expect an even more rapid dependence of B^ jB% 
near B2 = 0. As a consequence, it will be difficult to obtain the universal ratio 
B^ /B^ ~ 11.94 in calculations with finite range potentials. 

We now compare our results with the large- iV prediction Bn/Bn_\ 8.567 

from Ref. 0- Using the three-body results mentioned above, we find B^ / 'S3 = 
11.94. This number is considerably closer to the asymptotic value for Stv/Sjv-i 
than the value for N = 3: B^ / B2 = 16.52. The exact few-body results for 
N = 3, 4 are compared to the asymptotic prediction for S^t/Sat-i indicated 
by the dot-dashed line in Fig. ^ The dashed line gives an estimate of how the 
large- iV value should be approached. This estimate assumes the expansion: 



B N =(co + ^ + ^ + ..) 8.567^, (13) 
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leading to 



Bn 



8.567 + 0(N~ 2 ) . 



(14) 



B N -i 



The dashed line in Fig.^ was obtained by fitting the coefficients of the 1/N and 
1/N 2 terms in Eq. (|14|) to the data points for N = 3,4. A scenario where the 
corrections to Eq. 141) already start at 0{N~ l ) is disfavored by the data points 
in Fig. [I] 

Summary & Conclusion. — We have computed the four-body binding energies 
of weakly attractive bosons with a shallow two-body bound state in 2D. We 
found exactly two bound states and obtained the values B^ = 197.3(1) £>2 for 
the ground state and B^ = 25.5(1) B2 for the excited state. The four-body 
system is considerably closer to the universal large- N result B^/B^-i = 8.567 
derived in Ref. jS] than the three-body system. 

It would be important to test this prediction at even larger N both the- 
oretically and experimentally. Exact solutions of the quantum-mechanical N- 
body equations, however, are currently not feasible for N > 5. It would be 
interesting to see whether Monte Carlo methods such as the combined Monte 
Carlo/hyperspherical approach used for 3D 4 He clusters in Ref. ^5] can be 
adapted to this problem. This will not be an easy task, since the calculations for 
generic finite range potentials |12U13llTl| are not close to the universal zero-range 
limit. Finally, it would be valuable to create bosonic iV-body clusters in exper- 
iments and measure their size and binding energy. The possibility of creating 
such clusters in anisotropic atom traps should be investigated in the future. 
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